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A. Measurement Assurance Program for spectrophotometry is being established in order to assist laboratories 
involved in spectrophotoraetric calibrations. This paper deals with the preparation and calibration of neutral 
density glass filters for checking the linearity of photometric response, as applied to spectral transmittance 
measurements. Several sets of filters were prepared from suitable neutral glass to provide nominal transmit- 
tances of 92, 70, 50, 25, 10, 1, and 0.1% at a wavelength of 54S.5 nm. These filter sets will be available 
in three sizes: these are, 3& mm diameter aperture in 51 x 51 mm holder, 25 mm diameter aperture in 
51 x 28 mm holder, and 30 x 8 mm aperture in a cuvette holder. The filters were calibrated for spectral 
transmittance on the NBS Reference Spectrophotometer for high accuracy transmittanee measurements. 
Measurements were made with a 1.5 nm passband collimated sample beam. The filters were checked for 
uniformity and measurements were made to determine the effects of sample beam polarization. The transmil- 
tance data for the wavelength range of interest were analyzed by statistical methods to determine the ef- 
fects of passband for a range of 1 .5 nm to 10.5 nm passband. The results of these measurements are presented 
in tabular and graphical detail for the master filter set. 

Key words: experimental design; filter uniformity; linearity testing; measurement assurance program; neutral 
density glass; passband effects; polarization effects; polynomial fitting; spectrophotometry; statistics; transmit- 
tance standards. 



1. Introduction 

This paper deals with the preparation and calibration 
of neutral density glass filters for checking the linearity 
of photometric response, as applied to spectral transmit- 
tance measurements. 

The application of absolute techniques, such as the 
double-aperture method [l], 1 for checking the linearity 
of a spectrophotometer's response, is difficult to realize 
on many instruments because of problems involving sam- 
ple and reference beam geometry or detector configura- 
tions. The other approach to checking linearity is through 
the use of a set of transmittance standards. Several sets 
of glass filters have been prepared and calibrated for this 
purpose. A Measurement Assurance Program (MAP) is 
being established through which these sets of transmit- 
tance standards will be used for the purpose of improv- 
ing the accuracy of spectrophotometer measurements in 
laboratories participating in the program. 



"This project is supported by the Office of Measurement Services (Dr. B. C 
Belanger, Chief) and coordinated by Dr. L. J. Kieffer. 

""Center for Radiation Research, National Measurement Laboratory. 
tCenter for Applied Mathematics, National Engineering Laboratory. 
'Figures in brackets indicate literature references at the end of this paper. 



2. Selection ond Preparation 
of the Standards 

Some of the factors to be considered in selecting suitable 
filters for transmittance standards are: available range of 
transmittances, light scattering properties, uniformity, 
stability, passband sensitivity, and wavelength sensitiv- 
ity. Schott 2 neutral density glass types NG-4, NG-9, and 
NG-11 were selected on the basis of these characteristics 
and properties as well as previous experiences with Schott 
glass by NBS in the preparation of Standard Reference 
Material (SRM) 930 [2].* Two levels of transmittance 
were prepared from each of the glass types as follows: 



'Certain commercial materials are identified in this paper in older to adequately 
specify the experimental procedure. Such identification does not imply recommen- 
dation or endorsement by the National Bureau of Standards nor does it imply 
that the materials identified are necessarily the best available for the purpose. 

3 SRM 930 consists of three filters with nominal transmittances of 10%, 20%, 
and 30% calibrated at several wavelengths. Originally these SRM filters were 
designed for use by clinical chemists and are sold by the set. Our goal with the 
present filters was to provide a larger dynamic range so that the filters could be 
used to measure linearity in applications when other methods are difficult if not 
impossible. Also, these filters are not for sale, but are to be retained by NBS and 
issued only temporarily as part of a MAP measurement service. 
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Equivalent 
Glass Type Nominal Transmittance Transmission Density 



NG-11 


70% 


0.15 


NG-11 


50% 


0.3 


NG-4 


25% 


0.6 


NG-1 


10% 


1.0 


NG-9 


1% 


2.0 


NG-9 


0.1% 


3.0 



An additional filter of borosilicate crown glass having a 
nominal transmittance of 92% (transmission density 
0.036) is included in the set. 

The filter sets were prepared in three sizes in order to 
accommodate different instruments. The largest filters are 
mounted in 51 x 51 mm holders with a filter aperture of 
38 mm diameter. The intermediate size filters are mounted 
in 51 x 38 mm holders with a filter aperture of 25 mm. 
The smallest filters are mounted in cuvette holders and 
provide a filter aperture of 30 x 8 mm. 

The preparation of the filters from the stock glass was 
done in the NBS Optical Shop. The tolerances on 
parallelism of the two faces of the filters was maintained 
at approximately 0.01°. They were polished to a flatness 
of 3 fringes or better. 

Calculations of the required thicknesses were made 
from internal transmittance data supplied by the 
manufacturer. The internal transmittance of 1 mm path 
length for the three glass types is as follows: 



Glass Type 1 mm path 



NG-4 
NG-9 
NG-11 



0.32 

0.035 

0.78 



The values of transmittance for a wavelength of 500 nm 
were computed from these internal transmittance values. 
These transmittance values were converted to equivalent 
transmission density. Since the transmission density is 
linear with thickness, it is possible to determine, by 
graphical techniques, the required thickness for a 
desired transmittance value. Figure 1 shows the relation- 
ship of transmission density to thickness for the glass 
types NG-4, NG-9, and NG-1 1. The degree of accuracy 
in determining the required thickness for a given 



transmittance value by this graphical technique is 
dependent on the initial accuracy of the internal 
transmittance values used to determine the transmit- 
tance of 1 mm and 2 mm thicknesses. The actual 
transmittance values achieved by this technique closely 
approximated the desired nominal transmittance values 
for the filter set. 
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Figure 1. 



TRANSMISSION (%) 



Transmission density versus thickness of Sehott NG-4, 
NG-9, and NG-11 neutral glasses as derived Irom internal 
transmission data for 1 mm and 2 mm thicknesses. 



3. Transmission Analysis 

3.1 Light Scattering 

Scattering of a collimated sample beam transmitted 
through the neutral filters was quantitatively assessed by 
the ASTM recommended method [3] for measuring haze 
percentage. The measurements were made on the 
reference hazemeter [4] using the ASTM recommended 
geometry and methods. This instrument has a well col- 
limated circular incident beam having a color 
temperature of approximately 6800 K. A visual response 
filter at the detector modifies its response to give the 
hazemeter a peak sensitivity at approximately 550 nm. 
The instrument measures total forward scattering 
relative to the total transmission within the definition of 
this response and for the ASTM recommended 
geometry. The haze percentage determined for a set of 
filters was found to be as follows: 
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Filter 


Haze Percentage 


Nominal 
Transmittance 


1-1 


0.05% 


92% 


1-2 


0.02% 


10% 


1-3 


0.05% 


50% 


1-4 


0.06% 


25% 


1-5 


0.04% 


10% 


1-6 


* 


1% 


1-7 


* 


0.1% 



* The 1 and 0.1% transmitting filters were too low in transmittance 
for analysis by the hazemeter. It was concluded from these results that 
the filters scattered less than 0.1% of the transmitted sample beam. 



3.2 Survey Transmittance Measurements 

The purpose of a survey study of the neutral filter 
transmittance was to assess their over-all spectral prop- 
erties in the visible spectrum and select a specific 
wavelength for detailed analysis and calibration. Since 
the filters are not perfectly neutral, it was desirable that 
a wavelength be selected for calibration, at which the 
transmittance values of the filters would be least sen- 
sitive to variation in passband and wavelength scale 
errors associated with various spectrophotometers. 

With the exception of the borosilicate crown glass 
filter, the filters in the set are of a similar type glass with 
a common dye In three different concentrations for the 
NG-4, NG-9, and NG-11 glass types. The spectral 
transmission density curves of these filters in the spectral 
range between approximately 440 nm and 620 nm shows 
that the absorptions appear at the same wavelengths for 
the different concentrations associated with these glass 
types. Figures 2 through 8 illustrate the results of a 1-nm 
interval survey of the spectral transmittance of the filters 
between 440 and 620 nm. These measurements were 
made with a commercial spectrophotometer. The slit- 
width was maintained at approximately 0.1 mm (0.37 
nm passband). Further measurements were made of the 
10% filter with slit-widths of 1.0, 2.0, and 3.0 mm. 
These slit-widths correspond to passbands of 3.7, 7.4, 
and 11.1 nm, respectively. 

There are four wavelengths within the wavelength 
range 440 to 620 nm corresponding to peaks and valleys 
that could be used for the purpose of establishing a 
calibrated set of photometric scale standards. However, 
only one wavelength is needed and the absorption peak 
at approximately 548.5 nm was selected as being the 
most suitable. The spectral peak at approximately 464 
nm and the absorption valleys at approximately 510 and 
591 nm do not show a significant change in transmit- 
tance for the range of passbands used in this survey. A 



photometric resolution of approximately ±0.01% is re- 
quired in order to detect the effect of passband on the 
spectral transmittance measured at these wavelengths. 
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FIGURE 2. Spectral transmittance of Master Filter 1-1. 
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FIGURE 3. Spectral transmittance of Master Filter 1-2. 
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FIGURE 4. Spectral transmittance of Master Filter 1-3. 
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FIGURE 5. Spectral transmittance of Master Filter 1-4. 
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FIGURE 6. Spectral tranamittance of Master Filter 1-5. 
.0l05r 



650 




.0080 



400 



450 500 550 600 
WAVELENGTH (nm) 



650 



.UUIIU 










.00105 


- 




/ \ FILTER 1-7 


- 


iii 










Z .00100 

i- 

£ 

§ .00095 










r- 










.00090 








I 


.00085 




1 


i I 


1 



400 



450 



500 550 
WAVELENGTH (nm) 



600 



650 



FIGURE 7. Spectral transmittance of Master Filter 1-6. 



FIGURE 8. Spectral Transmittance of Maater Filter 1-7. 

3.3 High Accuracy Transmittance Measurements 

Calibration of the spectral transmittance of the master 
set of filters and three other sets was done on the NBS 
Reference Spectrophotometer [5,6,7] for measuring high 
accuracy transmittance. The measurements were made 
with a passband of 1.5 nm and a beam diameter of 10 
mm. Measurements were made for perpendicularly (s) 
and horizontally (p) polarized incident sample beams at 
1.5 nm intervals from 545.5 to 557.5 nm. Each filter was 
scanned for uniformity at 547 nm. The uniformity 
measurements were made at three locations; center, 2 
mm horizontally from center, and 2 mm below center. 
Transmittances for master set #1 filters are listed in table 
1. The results of the uniformity scan at 547 nm are listed 
in table 2. The transmittance values listed in table 1 are 
an average of the two polarizations. The values for both 
polarizations are listed for the uniformity measurements 
made at 547 nm. The instrumental uncertainties in the 
values listed in tables 1 and 2 are estimated to be 
±0.04% for filter 1-1 to ±0.0005% for filter 1-7. A 
complete description of the measurement sequence and 
data reduction for the high accuracy measurements and 
an analysis of the errors associated with such 
measurements is given in refs. [5-7]. 



4. Statistical Analysis 

A statistical analysis of the data obtained on the spec- 
tral transmittance of the neutral filters was used to deter- 
mine the magnitude of change in transmittance values 
due to changes in passband and also the magnitude of 
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TABLE 1. Spectral Transmittance of Master Set No. I for a 
Passband of 1.5 nm. 



Wave- 
















length 


Filter 


Filter 


Filter 


Filter 


Filter 


Filter 


Filter 


(nm) 


1-1 


1-2 


1-3 


1-4 


1-5 


1-6 


1-7 


536.5 


0.9164 


0.6902 


0.5177 


0.2339 


0.0955 


0.00919 


0.000915 


538.0 


.9165 


.6905 


.5181 


.2345 


.0959 


.00922 


.000922 


539.5 


.9165 


.6908 


.5185 


.2350 


,0962 


.00926 


.000925 


541.0 


.9165 


.6910 


.5187 


.2353 


.0964 


.00928 


.000929 


542.5 


.9166 


.6912 


.5190 


.2357 


.0966 


.00931 


.000933 


544.0 


.9166 


.6913 


.5192 


.2360 


.0969 


.00932 


.000935 


545.5 


.9166 


.6913 


.5192 


.2362 


.0970 


.00933 


.000936 


547.0 


.9166 


.6914 


.5194 


.2363 


.0971 


.00934 


,000938 


548.5 


.9168 


.6913 


.5193 


.2363 


.0971 


.00934 


.000938 


550.0 


.9167 


.6912 


.5191 


.2363 


.0971 


.00934 


.000940 


551,5 


.9167 


.6912 


.5190 


.2362 


.0971 


.00933 


.000937 


553.0 


.9167 


.6911 


.5190 


.2361 


.0970 


.00931 


.000934 


554.5 


.9166 


.6908 


.5187 


.2358 


.0968 


.00930 


.000934 


556.0 


.9168 


.6908 


.5184 


.2355 


.0966 


,00928 


.000929 


557.5 


.9168 


.6904 


.5180 


.2351 


.0964 


.00924 


.000923 



Table 2. Transmittance Uniformity of Master Set No. 1 at 547 nm. 

Location, Filter Filter Filter Filter Filter Filter Filter 

Polarization 1-1 1-2 1-3 1^ 1-5 1-6 1-7 

1, p 0.9166 0.6913 0.5193 0.2362 0.0971 0.00934 0.000939 

1, a .9166 .6913 .5193 .2362 .0971 .00935 ,000940 

2, p .9166 ,6914 .5193 .2361 .0971 .00934 .000937 

2, p .9167 .6914 .5193 .2362 ,0971 .00934 .000938 

3, p .9164 .6913 .5192 .2362 .0971 .00934 .000938 
3, p .9165 .6913 .5192 .2362 .0971 ,00933 .000939 



change due to location of the sample beam (filter 
uniformity). 

Preliminary data from the commercial instrument in- 
dicated that the central wavelength of interest was at ap- 
proximately 547 nm. However, a further analysis of the 
data from the high accuracy instrument indicated that 
548.5 nm was a better choice. 

Data are presented in table 3 to illustrate the effect of 
passband on the transmittance at the central wavelength 
for the master set #1 filters. These values are based on 
calculation for a triangular passband, using continuous 
integration. Figures 9 through 15 show these results plot- 
ted with the standard error limits. The data listed in 
table 4 illustrate the magnitude of change due to position 
location for the master set of filters. 

Remeasurement of the spectral transmiitances of the 
neutral filters was made one year after initial calibra- 
tions. The results of this second calibration for master 
set #1 are compared with the original calibration and 
listed in table 5 for wavelength 548.5 nm. Differences in 
spectral transmittances of the filters shown in table 5 are 
too small to be clearly interpreted as changes with the 
possible exception of filter 1-1. Here the apparent 
change in transmittance was -0.0009. Further 
measurements at longer time intervals will be required to 
confirm any real changes in spectral transmittance due 
to aging. 

Additional details of the statistical analysis are given 
in appendix A. 



Table 3. Transmittance of Master Set No. 1 at 548.5 nm as a Function of Passband for a Triangular Passband. 



Passband 
(nm) 



Filter 
1-1 



Filter 
1-2 



Filter 
1-3 



Filter 
1-4 



Fitter 
1-5 



Note: uncertainties (one standard deviation) are in parentheses 



Filter 
1-6 



Filter 
1-7 



1.5 


0.91665 

( 31 


0.69133 

< 2) 


0.51928 

( 2) 


0.23632 

( 2) 


0.09713 
{ 1) 


0.009342 
( 2) 


0.000939 

( 4) 


3.0 


.91665 
( 31 


.69132 
( 2) 


.51927 
( 2) 


.23632 

( 2) 


.09711 
( 1) 


.009341 

( 2) 


.000939 

( 4) 


4.5 


.91665 

( 3) 


.69130 

( 2) 


.51924 
( 2} 


.23628 

( 2) 


.9709 
( 1) 


.009338 

( 2) 


.000938 

( 4) 


6.0 


.91665 
( 3) 


.69128 

( 2) 


.51921 

( 2) 


.23624 
1 2) 


.09706 
( 1) 


.009335 

{ 2) 


.000938 

f 4) 


7.5 


.91665 
( 3) 


.69125 

( 2) 


.51916 

( 2) 


.23618 

( 21 


.09703 

( 11 


.009332 

( 2) 


.000937 

( 4) 


9.0 


.91665 

( 3) 


.69121 

( 21 


.51911 

( 21 


.23616 
( 2) 


.09698 

( 1) 


.009327 
( 2) 


.000937 
( 4) 


10.5 


.91665 
( 3) 


.69117 

( 2) 


.51905 

( 2) 


.23604 

< 2) 


.09693 

( 1) 


.009321 
( 2) 


.000936 

( 4) 
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FIGURE 9. Spectral transmittance at 548.5 nm versus bandpass for 

Master Filter 1-1. 
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FIGURE 12. Spectral transmittance at 548.5 nm versus bandpass tor 
Master Filter 1-4. 
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FIGURE 10. Spectral transmittance at 548.5 nm versus bandpass for 
Master Filter 1-2. 
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FIGURE 11. Spectral transmittance at 548.5 nm versus bandpass for 
Master Filter 1-3. 
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FIGURE 13. Spectral transmittance at 548.5 nm versus bandpass for 
Master Filter 1-5. 
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FIGURE 14. Spectral transmittance at 548.5 nm versus bandpass for 
Master Filter 1-6. 
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Figure 
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15. Spectral transmittance at 548.5 nm versus bandpass for 
Master Filter 1-7. 



TABLE 4. Magnitude of Filter Nonuniformity 



Filter 
No. 



1-1 
1-1 

1-1 



Location 



Transmittance 

at 548.5 nm 



Location 
Difference 



0.91658 U)-(2) = -0.00008 

.91666 (2M3) = .00018 

.91648 <1M3) = -00010 



92%. These filters will be used in a Measurement 
Assurance Program (MAP). The purpose of this pro- 
gram will be to assist laboratories wishing to maintain a 
high level of confidence in the accuracy of their spec- 
trophotometric measurements. 

Data are presented that show the spectral transmit- 
tances of these filters at 548.5 nm to be relatively insen- 
sitive to variations in passband and slight errors in in- 
strument wavelength setting. The filters are of sufficient 
uniformity for use as transmittance standards. 



Table 5. Spectral Transmittance of Master Set No. 1 for a 
Passband of 1.5 nm at 548.5 nm. 





(Apparent changes in transmittance) 






First Calibration 


Second Calibration 


Apparent 


Filter No. 


(Oct. 1980) 


(Oct. 1981) 


change 


1-1 


0.9168 


0.9159 


-0.0009 


1-2 


.6913 


.6914 


.0001 


1-3 


.5193 


.5196 


.0003 


1-4 


.2363 


.2365 


.0002 


1-5 


.09712 


.09722 


.00010 


1-6 


.009337 


.009350 


.000013 


1-7 


.0009376 


.0009395 


.0000019 



1-2 
1-2 
1-2 



.69132 
.69138 
.69131 



U)-<2) = 
<2)-(3) = 
<l)-<3> = 



.00006 
.00007 
.00001 



6. References 



1-3 
1-3 
1-3 

1-4 
1-4 
1-4 

1-5 
1-5 
1-5 

1-6 

1-6 
1-6 

1-7 

1-7 
1-7 



1 

2 
3 

1 
2 
3 

1 

2 
3 

1 
2 
3 

1 

2 
3 



.51926 
.51928 
.51919 

.23624 
.23615 
.23619 

.09709 
.09708 
.09708 

.009341 
.009342 
.009331 

.000939 
.000938 
.000937 



UM2) 
(2M3) 
UM3) 

(11-121 
<2)-(3) 
(l)-(3) 

(l)-(2) 
(2)-(3) 
(l)-(3) 

UM2) 
(2M3) 
(1H3) 

(l)-(2) 
(2M3) 
UM3) 



.00002 
.00009 
.00007 

.00009 
.00004 
.00005 

.00001 
.00000 
.00001 

.000001 
.000011 
.000010 

.000001 
.000001 
.000002 



5. Summary 

One master set and three working sets of neutral den- 
sity glass filters have been calibrated for use as spectral 
transmittance standards for checking the photometric 
scale linearity of spectrophotometers. Each set consists 
of seven filters ranging in transmittance from 0.1 to 
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Appendix 

Statistical Design and Analysis for Calibrating the 

Transmittance of Filters for Spectrophotometer 

Linearity Testing 

For four sets of filters, each set containing seven filter 
types of various transmittances, an assigned transmit- 
tance value needed to be determined when a triangular 
passband of given width is superimposed. Additionally, 
the polarization and location of the spot size on the filter 
had an unknown influence on the transmittance value. 
Although a commercial instrument can scan the entire 
spectrum at little cost, the measurements are an order of 
magnitude less accurate than those obtained using the 
high accuracy instrument. For this reason, a design was 
developed to determine: 

(1) The values to be measured for purposes of 
calibration, and 

(2) The existence of a possible effect due to polarity 
or the effect of location on the filter. 

The motivation behind the choice of design is discussed 
in section A below. In the event of (2), the calibration 
must be reconsidered. On the basis of these extremely 
precise measurements, an effect due to location was 
statistically significant. This suggests that either: 

(1) The calibrated values can be used only for the 
central location on each filter; or 

(2} Additional measurements will be required to 
determine the magnitude of changes in transmit- 
tance due to spot location. 

Since the observed differences are small «0.2%; see 
table 4), the location effect may not interfere with prac- 
tical usage of these filters. However, measurement as 
near to the central location on each filter as possible is 
recommended. The conclusions regarding location and 
polarity effects are presented in section B, and the 
estimation of the transmittance curve leading to the 
calibration values is shown in sections C and D. The 
method of extension for general passbands is given in 
section E. 

A. The Design 

Most organizations using these filters will have 
triangular passbands on their measuring equipment. 
Therefore, adequate determination of the transmittance 
using this passband is essential. However, other pass- 
band shapes may be used, and for these, a more general 



approach is taken. The problem then becomes one of 
estimating the transmittance curve with high precision. 
Since the measurements on the high accuracy instrument 
are time-consuming, a design was chosen so as to permit 
adequate curve estimation with relatively few points. 

1. Preliminary analysis: Choice of Wavelength. 

The choice for a wavelength to be calibrated required 
consideration of: 

(i) Varying passband widths, and 
(ii) The potential for imprecise wavelength specifica- 
tion. 

We consider neighborhoods of 21 nanometers, since the 
widest passband that is most commonly used weights 
frequencies within a range of this length. The chosen 
wavelength is: 

A„ = 548.5 mn. 

This subsection explains this choice in light of considera- 
tions (i) and (ii). 

For a passband having width w and functional form 
/i„(A), the transmittance at wavelength A is given by the 
convolution integral 

t„U )=/ UA)/(A-A o )dA / J UA)dA, 

where /(A) is the spectrum (transmittance curve). If the 
passband is normalized so that 

A„+Q) 

/ hJX)dX = 1, 

A„-o> 



then this simplifies to 



t„(A ) = J hJX\f(X-XJdX . 

A-ta 



(1) 



In all that follows, we will assume that the passband ft w 
has been normalized in this way. 

In order for this average value to be relatively insen- 
sitive to wavelength specification, /(A) needs to be a 
smooth, slowly-varying function. In mathematical 
terms, we search for a neighborhood where the first two 
derivatives of the spectrum are fairly small. This sug- 
gests a region where /(A) is nearly constant, or, at worst, 
a quadratic having very small degree of curvature. A 
locally linear spectrum with a noticeably nonlinear slope 
would yield a seriously biased estimate, whereas a locally 
quadratic spectrum, centered at a peak or trough, would 
have a negligible linear term and therefore a bias which 
depends primarily on the second-order term {hopefully 
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small). Preliminary readouts on a commercial instru- 
ment suggested a central wavelength of 547.0 nm for 
most filters. On the basis of the first set of measurements 
on the high accuracy instrument, however, this was later 
amended to A„ = 548.5 nm. 



3. Location of Points 

We consider two criteria in selecting the four 
wavelengths (two symmetrically placed on either side of 
the central wavelength): 



2. Number of Points 

Assuming that the spectrum is locally quadratic and 
has a maximum around A = 548.5 nm, our goal is to 
estimate the three parameters a, ji, y in the approxima- 
tion: 



fiX) = a + jiX + Y k 2 - 



(2) 



In order to estimate a, /?, and y, by a, j5, y, the proposed 
design specifies a total of five points, one at the central 
wavelength, and one each at equal distances ±26 from 
the central wavelength. 

Let x — A - 548.5, so as to center the relation at the 
origin. A minimum of three points is required to specify 
the parameters. This leaves no room for the assessment 
of error. With more points, the variance of (1) will be a 
linear combination of the variances and covariances of 
the parameters, and will be dominated primarily by the 
variance of the constant term. If our chosen points are 
symmetrically placed about 0, it can be shown that: 



where 



Var (a) =oV(n-cVe), 



Xi >0 
X t X3 



(3) 



As (Z*, 4 ) < CZxfY [x i > 0), the smallest value of cVe is 2 
{all points at the origin except for the two at ±1). A 
lower bound on (3) is therefore: 



Var (5) = oV(u-2). 



(3*) 



If we concern ourselves with only the first term, then, 
relative to o 2 , eq (3 ') gives the reduction in the overall 
variance that we may hope to gain in our variance of the 
transmittance value. With n = 3, there is no reduction; 
with n = 5, the variance is already reduced to 33% of its 
value. An extra two points provides only an extra 13% 
improvement. Thus, the design specifies a minimum of 
two points on either side of the central wavelength. 



(i) Small error in quadratic interpolation of the 

transmittance curve; 
(ii) Small variance in the transmittance estimate 

given by (I). 

On the basis of these two criteria, the proposed design 
specifies measurements to be taken at: 



A„, K 



± 6 nm, X a ± 9 nm. 



The reasoning behind this choice is explained by con- 
sidering (i) and (ii) separately. 

(i) Error is quadratic interpolation. A bound for the 
estimate of error in polynominal interpolation of degree 
n at the points a < x u x 2 , ... x n+l < b is 



|e(A)| < max |f< n+1 >U>| 
a<A<b 



|A-*i[ X~x 2 \ ■ • ■ \X-x„, 



(see, e.g., [8]). The first part of this error depends upon 
the transmittance curve and led to the choice of J a = 
548.5 in subsection 1 above. Now we concentrate on 
selecting x x , ... , x a+1 so that the error is as small as possi- 
ble. The answer to this problem is given by the zeroes of 
the Chebyshev polynomials, namely 



n + 1. 

(4) 



** = <10.5nm)'cos [(2ft-l>/(2n+2)], ft = 1, 



For the problem at hand, n+1 = 5; thus: 

xi = -X5 = 9.986 nm from A 

Xi = -Xi =6.172 nm from A, 

*3 = nm from A„ 



As most passbands are in multiples of 1.5 nm, the closest 
multiples to these points are ±6 nm and either ±9 nm or 
±10.5 nm. However, the nature of the triangular pass- 
band which is frequently used in practice would assign 
zero weight to the values ±10.5 nm. Since a primary 
goal is the estimation of the transmittance using this 
passband, we choose ±9 nm. 
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(ii) Minimum variance of transmittance estimate. If 
we estimate the transmittance curve using a quadratic 
function, viz. 

/<*)=<? +Px+yx?, <5) 

then the transmittance using passband of width w is 

f w (A ) = f hJw+XJ i& +px + yx 2 )dx. 

Let 

c,= f xihjx) dx 

- uj 

Then 

tJA„) =ac + /3c t + yc 2 = c'0 
where 

c' = (c„, Ci, c 2 ) 

0' = 6,/3,y). 

Then minimizing the variance of c '0 is equivalent to 
min Var(c'0) = min c' Var(0) c. 

A design which minimizes the variance-covariance 
matrix Var(0) of the parameter estimates is given in 
[9], The design would place the two observations each at 
A + 10.5 nm in addition to the one at A . If our function 
is truly quadratic, such a design would be optimal. 

For many reasons, however, we modify the optimal 
design which permits more flexibility in our choice of 
model. Such a design is suggested in [10]. The design 
recommends two different values rather than repeating 
them at the endpoints. 

Closely related to this design is the one which 
minimizes the maximum variance of the best linear un- 
biased estimate of the function f(x) given by eq (5). 
Reference [11] shows that the five points should be 
placed at the zeros of the polynomial 

10.5 (1 - x 2 } PAx) , 

where Pi ' is the derivative of the fourth degree Legen- 
dre polynomial 

P i '(x) = n.5x i -1.5x. 



Hence, the five points are 0, ±7,1, ±10.5 nm from A„. 
Again, since our passband applies decreasing weight to 
f(x) as x is further from the origin, we choose to make 
these points in towards 0, to ±6 and ±9 nm. 

4, Tests for Polarity and Spot Location 

Additional measurements at the central wavelength 
are needed to provide tests for difference in polarization 
and spot location. In addition to the measurement at the 
central wavelength with polarity 1 (point 1 at location 
(0,0)), measurements were taken at 

Point 2: Spot location (2,0), Polarity 1 
Point 3: Spot location (0,2), Polarity 2 
Point 4: Spot location (0,2), Polarity 1. 

This permits a check for a difference due to polarity by 
comparing the third and fourth points, and check for 
location differences by comparing points 1 and 4 (or 1 
and 2 or 2 and 4). A more extensive check based on a full 
2X3 factorial combination was made on the master set 
#1. 

B. Results: Polarity and Location 

The reported results are based on measurements taken 
on four sets of filters. A complete series of 15 
measurements across the 21-nm range at 1.5 nm spacing 
was taken on the master set, as a check for the adequacy 
of the five-point design used on the other three sets. In 
addition, polarity and location were tested on each filter 
in all sets. The results of these tests are reported in this 
section; transmittance calculations are reported in sec- 
tion C. 

For set 1, six measurements were taken at A = 547.0, 
at both polarity 1 and 2 at the center of the filter (0,0) 
and at placements of two units to the right (2,0) and 
above (0,2) the center. This permits a 3 X 2 factorial 
analysis for detecting differences due to location and 
polarity. 

For sets 2, 3, and 4, Mests on location and polarity 
were calculated. The test on location was deemed signifi- 
cant if the comparison of either (0,0) and (0,2) or (0,2) 
and (2,0) indicated a chance of less than 1% under the 
hypothesis of no effect. 

None of the tests of polarity was significant at the 0.01 
level. (One filter from one set gave significance of 0.025. 
Out of 28 such tests, the chance of obtaining one or more 
spurious significances is more than 50%, so this is hardly 
surprising.) For location, however, the following sets 
showed a significant difference (level given in 
parentheses): 



34 



Filter 


Set (level of significance)* 


2 


3. (0:0004* 


3. 


2 (O'.O'ISK 3- (O'.O02>, 4K0>.0OM)i 


4 


2' (O.0M3)i,. S (O'.OW, 4 (OiM> 


5 


3 (O.O'IK 4 (O-.O0O3)- 


6- 


1 (0.0061, 4 (O.OOM). [3 (O-.Oa-)ll 


7 


1 (0.0002k 4 (&.002S)' 



Out of 28 tests of signffieaace, the- chance of 131 or more 
coming up* significant at the Ov0»I level when in fact loca- 
tion! is irrelevant is; 



2» 
Jfc=13i 



l 2 h tMlt L99F* =0 



and' the chance that at least 2! of the 4 sets would show 
significance on a given filter type- is. 



4 

2 

ft=2 



(j* j (o. 01)' (0i.99f"* = 0.0006. 



Consider one filter type from each set having a 
nominal transmittan.ee teay, 69%.). The average value of 
Iransmittance, given by the parameter a, is likely to be 
specific to each filter in the different sets. However, the 
curvature parameters, fi and y are likely to be common to 
all 4 filters, having a given nominal transmittance. Thus 
we use all these observations, to fit a relation of the form 

f{ XlJ \ = otj. + p-xu + yxl , (7) 

i ] — 1, ..,, tij {= # of measurements in jth set). 
/ = 1, ..., 4{= # of sets). 

Note that 

n, = 21, 
{x„\ = {-12., -10.5 7.5, 9, -1.5 -1.5} , 

since set 1 has the full 15 measurements and 6 additional 
measurements for polarity and location, and 

nj = 8, {*„} = {-9, -6, 0, 6, 9, 0, 0, 0} 



For 7 types of filters, the overall 1 level 1 of significance is 
approximately 0.004* 

These tests of significance suggest that it is highly 
unlikely that location on the filter in measuring transmit- 
tance is irrelevant. However, all tests were based on the 
internal standard errors of the four measurements made 
within the twenty-minute measurement period. This 
standard error reflects only the error of the four internal 
readings but does not reflect the measurement-to- 
measurement variation caused by apparatus set-up, 
filter placement in the wheel or polarity switches. As 
such, the actual error in taking successive measurements 
may in fact be larger than the reported standard error. 
However, measurements on these filters as near to the 
central location as possible is recommended. The 
transmittance values using triangular passbands have 
been calculated assuming location is irrelevant (thereby 
using all data in the estimation of the transmittance 
curve) and assuming location has a significant effect on 
the transmittance value (using only the data at the center 
(0,0) spot). Only the latter are given for the master Set #1 
in table 3 of the main report. 



C. Estimation of Transmittance Curve 

As there are four sets which contain all filter types, it 
will be helpful to use all filters having the same nominal 
transmittance in estimating the transmittance curve. We 
fit quadratic functions of the form in eq (5). 



for sets 2, 3 and 4. Note also that the average level, aj, is 
specific to the filter in the set, as the level may reflect the 
amount of dye that is contained in the filter. However, 
the parameters p and y are likely to be influenced by the 
properties of the dye contained in this glass and are thus 
common to the filters with the same nominal transmit- 
tance from all four sets. This joint estimation permits the 
more accurate calculation of {S and y in sets 2, 3, and 4, 
where the number of observations is only eight. 

Since transmittance may depend on location (section 
B), we also estimate a, (i, y, for each of the 7 filter types, 
using only those observations at location (0,0). In this 
case, 

ni = 17, {*„} = {-12, -10.5 . . ., 7.5, 9, -1.5, -1.5} 
nj = 5, {*„} = {-9, -6, 0, 6, 9}, / = 2, 3, 4. 

These parameters were not found to differ from those 
using all the data by more than two standard errors. In 
light of potential location differences, however, table Al 
provides these estimates of the parameters a, fi and y for 
the master set #1 filters. 

D. Estimation of Transmittance Values Using 
Triangular Passband 

As indicated in eq (1 }, when the transmittance curve is 
superimposed with a normalized passband centered at 
A , the observed transmittance can be expressed as a con- 
volution integral. Using a triangular passband of width 
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w, symmetric about X , the result is a weighted average 
of all transmittances around the centra] wavelength. 
This weighted average may be computed either as a 
discrete sum: 

T w (A a ) =X h„(\i)yi yi = tnmsnrittancc at wavelength A;; 

=2 h w (x<+K)y;, *i = i.- - A„, (8) 

where 

hjx) = uH[l - sign(;e)-.%/tti], -w < x < tc, <9) 

or as an integral as in eq (I). The first approach, based 
on eq (8), is a discrete computation based on the ob- 
served values y,. The second approach, eq (1), uses the 
values y, to estimate f(x) (eq (5)) and integrates directly. 
For sets 2, 3, and 4 where data inside ±6 nm exists only 
at the central wavelength, eq (1 ) is clearly the method of 
choice. Furthermore, the estimation of/fccj uses all of the 
data in its estimation (section C), and is therefore likely 
to remove much of the variability in the values y,-. 

The standard error of t„ can be computed using the 
reported standard errors <s ; ) associated with each 
measurement: 



SE ( tJ = nhJxW sf. 

Using the triangular passband of width w (eq (9)) and 
the fitted relation (7), we have that the transmittanee of 
a filter from the jth set is 



UK) = f 



hM ■ fcj + Ptt-h) +yU-A,)] 2 dX 



ccj + (tc 2 /6) ■ y 



(10) 



for which 



S£(tJ =[Var(S,) + 
= SE (aj); 



Table Al. Coefficients /or Estimation of Transmittanee Curve for 
Set If I (includes data jrom center location only). 



i(x) = w + j9x+fx 2 x = wavelength -548.3 
Table gives constant term; all standard errors given in parenthe 





a 


i? X 10 s 


yX 10 1 


Filter 1 


.9166325 


1.6493 


-0.002562 




( 200) 


( .1877) 


■( .032669) 


2 


.6913352 


-2.0067 


-0.91372 




( 117) 


(1.50301 


( .02616) 


3 


.5192856 


-2.7205 


-1.2973 




( 135) 


( .1151) 


( .02001 


4 


-2363304 


H-l-1966 


-L5858 




( 147) 


( .0923) 


! .0161) 


5 


.0971290 


+0.8657 


-1.0699 




( 068) 


( .0736) 


( .-0128) 


6 


.00934226 


-.05959 


^0.11355 




( 113) 


(.00979) 


( -00170) 


7 


.000939014 


-0.008245 


-a 01 7568 




( 266) 


< .002893) 


( .000504) 



[w;V36)Var<y) + <u;V3)Cov(S,,y)]y 2 



since the variance of the eurvature parameter y and the 
covariance term are typically two orders of magnitude 
smaller than the variance of the fitted constant term. 

E. Other Passbands 

Table 3 of the main report gives transmittances for the 
master set #1 filters for a triangular passband. For 
passbands other than those listed in table 3, the 
transmittanee may be calculated directly via eq (10), us- 
ing the estimated coefficients eK, ft, and y listed in table 
Al. 
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